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Abstract 

We present an exactly solvable model for queue dynamics. Our model 
is very simple but provides the essential property for such dynamics. As 
\Q . an example, the model has the traveling cluster solution as well as the 

C> homogeneous flow solution. The model is the limiting case of Optimal 

■ Velocity (OV) model, which is proposed for the car following model to 

' | ■ induce traffic jam spontaneously. 

o. 

The concept of queue dynamics offers simple 1-dimensional models for socio- 
economic and complex multi-body physical systems, such as traffic and granular 
transport problems. Recently, several approaches have been developed on this 
field |JXL 0, ^ ]5], [/J. We propose an exactly solvable simplest model of this kind. 
Our model describes the general aspects of queue dynamics, and it can be widely 
applicable to such systems. But the following discussion is presented with the 
terminology of traffic problems. 
The model is 



x n = a { V(Ax n ) - x n } , (1) 



where 



A3^n %n— 1 *£n> (2) 

for each car number n (n = 1, 2, • • • ). x n is the position of the n-th car, Ax n is the 
headway of that car. Dot denotes the time derivative, a is a sensitivity constant, 
which we set the same value for all drivers. The function V(Ax n ), which is called 
"OV-function", is 

V(Ax n ) = v max ■ 9(Ax n - d), (3) 

where 9 is a Heaviside function. It decides the optimal velocity (the safety 
velocity) according to the headway: (a) if the headway is less than d, a car 
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should stop; (b) while the headway is larger than d, a car can accelerate to move 
with the maximum velocity v max . 

For both cases, the movement of each car is easily derived as follows, 
fa) In the case of Ax n < d 



x n (t)=x n (t ) + ^{l-e- a ^}, (4) 

with the initial condition, t = t , x(t ), x(to). 
(b) In the case of Ax n > d 

*»(«) = s B (f ) + " *o) - ^"^ {1 - e- a ^}, (5) 

a 

with the same initial condition. 

For the purpose of obtaining the solution of jam flow, we should consider the 
two basic processes of car moving. Suppose a jam exist in the lane, a car moves 
from "free driving region" into jam, and another car escapes from jam to free 
driving region. These two processes can be expressed by the above two solutions 
with appropriate connection conditions. We assume the ideal case: cars stop in 
jam with the same distance axj(< d), and move at the maximum velocity v max 
with the same distance AXp(> d) in free driving region. 

First, we investigate the process of car moving from jam to free driving region. 
When the headway of the top car in a jam cluster is d, we set t = to. At this 
time the position of the car is set as x = 0. The n-th car in the jam moves as 
the following formulas. 

x n {t) = -n ■ AXj (t < t < t n ) , (6) 

and 

x n (t) = -n-AXj + v max {(t-t n )-^-[l-e- a{t - tn) ]} {t n <t), (7) 

where t n is the time when the headway of the n-th car is d. After t n the car 
begins to move and escapes from jam. We note t < t 1 < ■ ■ ■ < £ n _i < t n < ■ ■ ■, 
and t n is defined by 

Ax n (t n ) = d . (8) 
From and (|7]), the definition (^) is written as 

axj + v max {{t n - t n - x ) - -[1 - e -°(*»-*»-0]} = d. (9) 

a 

We have derived the sequence of equations for the definitions of t n (n = 1, 2, • • •), 
which has the solution for t n > t n -i- It is easily obtained by putting r = t\ — 1 = 
■ ■ ■ = t n — £ n _i = •••(> 0), and that is the unique solution, which is given by the 
following equation 

ax j + v max {r - -(1 - e- aT )} = d . (10) 

(X 
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It is easily seen the velocity of the n-th car, x n (t) of (|7|), converges to v max 
for sufficient large time (t 3> t n ). This means the car reaches the free driving 
region, and the headway of this car becomes ax f , which is expressed as 

lim Ax n (t) = axf , (11) 

t — >oo 

with using (|7|). Thus, we have obtained the following simple relation. 

AXj + V max T = AX F . (12) 

Next, we make an analogous investigation for the process of car moving from 
free driving region into jam. When the headway of the car positioned just behind 
a jam cluster is d, we set t = to and the position of the car xq = 0. The n-th car 
in the free driving region moves as the following formulas. 

x n (t) = -n- Ax F + v max (t-t ) (t <t<t n ), (13) 

and 

x n (t) = -n-Ax F + v max {(t n -t ) + ^[l-e- a(t - t " ) ]} (t n <t), (14) 

where t n is the time when the headway of the n-th car is d. After t n the car begins 
to decelerate and moves into jam. We note t < t\ < • • • < t n -\ <£„<•••, and 
t n is defined by 

&x n (t n ) = d . (15) 
From (^) and ([14]), the definition (Tl5| ) is written as 

AXp + v max {-{t n - t ft _0 + -[1 - e -«<*»-*»-i)]} = d . (16) 

(X 

As the same as the previous case, we can obtain the unique solution of the 
sequence of the equations for t n with the conditions, r = t\—to = ■ ■ ■ = t n —t n -\ = 
■ • ■ (> 0), where r is given by 

Ax F + v max {-r + Ul-e- aT )} = d. (17) 



The velocity of the n-th. car, x n (t) of (|I4]) , converges to for sufficient large 
time (t ^> t n ), which means the car reaches a jam cluster and stops. So, the 
headway of this car becomes axj, which is expressed as 

lim Ax n (t) = axj , (18) 

t — >oo 

with using (|i~4|). Again, we have obtained the same relation as fll2D , which means 
the value of r is the same for both processes we have discussed above. 

Now we can solve the equations (fLOl), (|l|) and fllTD , thus ax f , axj and r can 
be expressed using a, d and v max as follows, 



3 



axf = d + 
ax j = d — 
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t is determined by the following equation 



ar = 2(1- e" 



(19) 
(20) 

(21) 



which has the solution ar ~ 1.59. As the result, each car is moving just the same 
manner as its front moving car with the time delay r, which is induced in the 
process of obtaining the jam flow solution. And it is proportional to the inverse 
of sensitivity 1/a. The collective movement of these cars forms a jam cluster. 

In Fig.|l] the moving of several successive cars are shown. Their orbits are 
consisted of (|13|), (0), © and (^) with the time delay r. Fig.|| shows the velocity 
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Figure 1: An example of jam flow. The movement of five successive cars 
are drawn by bold lines, which forms a jam cluster. The dashed lines show 
the movement of a jam cluster, which velocity is given by (22) or (23). 



of car with time development, which is easily obtained as the time derivatives 
of these formulas. Rigorously say, it takes infinite time for car velocity to reach 
Vmax or 0. Practically, an appropriate finite time is enough to be considered as 
infinite, as you see in Fig.|^. 

In Fig.|T] a jam cluster moves backward against the direction of car moving. 
The velocity of the cluster is defined by the moving of the front position of the 
cluster, which is obtained by (^ or (0) as 

Vjam = " — " = — • (^J 

"ti ^n— 1 ' 
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The same result is given by the moving of the rear point of the cluster from (|TB|) 
or flljD with fljjl). Using the values of axf and axj, the velocity is also defined 
as H 

AXp ■ V(AXj) — AXj ■ V^AXf) AXj 



J jam 



AXf — AXj 



(23) 



We obtain the same result as (|22|). We note each cluster should have the same 
velocity given above, if a jam flow contains several clusters. (Actually, you can 
see such a case later in Fig.|l.) 
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Figure 2: The changing of car velocity of three successive cars are drawn. 
The maximum velocity is set as v max = 2. We can see a kink like shape. 



The profile of a jam flow is clearly described as the trajectory in the phase 
space of headway and velocity (ax, x). In order to draw this, it is enough to 
check the movement of two successive cars. The movement is consisted of two 
basic processes: the moving from jam to free driving region and the opposite 
case. 

First, we check the process of car moving from jam to free driving region. 
The moving of the (n — l)-th and n-th cars is divided into three stages presented 
in Table |l|. 



Table 1: The movement of two successive cars from jam to free driving region. 



t 


to 


tn—l 


t n 




ax„_i < d : (|) 


Ax n _i > d : © 


•En 


AX n < d 


(i) 


Ax n >d : © 



i) to < t < t n -i : The headway of the n-th car is ax j and its velocity is 0. The 
n-th car stays the point (ax j,0) in the phase space in this period, which means 
the car stays in jam. 
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ii) t n -i <t<t n : The velocity x n is 0. The headway is 

Ax n = axj + v max {(t - * n _i) - -[1 - e -«(*-*»-i)]} . (24) 

(X 



From (P0|) and (|2lD , the headway changes axj < ax„ < d for i n _ 1 < t < t n . 

Thus, the trajectory is the line (axj, 0) - (d, 0). 

hi) t n < t < oo : In this period, the headway and velocity are 

Ax n = AXj + v max {(t n -t n ^ 1 ) + -[~e' a{t - t - ) + 6-^-^}} , (25) 

a 

x n = v max {l - e-^-^} . (26) 
From these equations with using (|20|) and (pT|) , we derive the following relation 

TX 

Ax n = d+ — (0<x n < v max ) . (27) 

The trajectory is the line (d, 0) - (axf, i> maa; ). 

Next, we turn to the car moving from free driving region into jam. The 
process is also divided into 3 stages, which is presented in Table §. 

Table 2: The movement of two successive cars from free driving region into jam. 



t 


to 


t n -i 


tn 






ax„_i > d : (|T5) 




AX n _i < d 


: (14) 




•En 


AX n > d 


am 




AX n < d : 


(M) 



i) t < t < t n -i '■ The headway of the n-th car is AXp and it's velocity is v ma x- 
The n-th car stays the point (AXF,v max ), which means the car moves in free 
driving region. 

ii) t n _i < t < t n : The velocity x n is v max . The headway is 

AX n = AXf + V max {-(t - t n -l) + - [1 - e^-^]} . (28) 

From (|H]) and (^Tj), the headway changes AXp > Ax n > d for t n _i < t < t n . 

Thus, the trajectory is the line (axf, v max ) - (d,v ma x)- 

hi) t n < t < oo : In this period, the headway and velocity are 

AX n = AX F + V ma A-(tn-t n -l)-^l-e~ a{t ~ t " ) +e- a(t - t "- l) ]} , (29) 
X n Umax^ ^ ^ ■ (30) 

From these equations with (^) and Q2"Tp, we derive the following relation 

TX 

AX n = AX j + — - (V m ax > X n > 0) . (31) 
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Figure 3: The hysteresis loop of jam flow together with the OV- function 
V(ax). Each car moves along this loop in the direction of arrow with time 
development. 



The trajectory is the line (d, v max ) - (axj,0). We summarize the above results 
to Fig.^. The car movement in jam flow solution is represented as the square 
shaped closed loop in the phase space. All cars are moving along this loop in 
stable, which can be understood as a limit cycle. 

In our previous work, we proposed the car following model (OV model) 0, 
whose OV function is hyperbolic tangent instead of (|3|). We found by the nu- 
merical simulation the model organized the limit cycle figure called "hysteresis 
loop" , which is attractive in the whole phase space (3|. The difference between 
our simple model and the previous case is only the shape of hysteresis loop. We 
are convinced that the essential property for the dynamics of jam flow solution 
of OV type model is clearly understood in our simple model. 

The size of hysteresis loop AXp—AXj = v max r, which determine the amplitude 
of jam cluster, is characterized by the induced time delay r. It is proportional to 
the inverse of sensitivity 1/a. The loop shrinks to the vertical line ax = d, (0 < 
x < v max ) as a — > oo, in which case jam does not appear. 

We should note the all above results do not depend on the total number 
of cars and the length of lane. We do not need to set the periodic boundary 
condition on the lane. If we set them as some conditions, it is easy to calculate 
the total number of cars in jam (the total sum of cluster size) or how long a car 
stays in jam. 

Finally, we check the above results against the simulation data. The aim is to 
know how well the analytic solution is realized, which we have obtained by the 
infinite time approximation. We set the parameters of the model as a = 1, d = 2 
and v max = 2. In this case axf — 3.59 and axj ~ 0.41 are derived from ( |T!3| ) and 
(|20|). The simulation is performed with putting cars on the circuit, whose length 
L = 200, and the total number of cars iV = 100. 

Fig.[| is the plot of the position of all cars with time development. The 
initial condition is set as all cars are uniformly distributed with the distance 
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ax = L/N = 2(= d) and move with the same velocity x = 1. In this case, the 
initial movement is highly unstable. We can observe the growth of jam clusters, 
and they are moving backward in stable with the same velocity, which value is 
in agreement with the analytic result (p2]). 
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Figure 4: The simulation data of the positions of all cars in the circuit 
with time development. The solid line shows the orbit of a sample car. 

Fig.[| is the snapshot of headway distribution of the jam flow on the circuit. In 
this case cars are gathered to one big stable cluster, which simulation is performed 
with a different initial condition from the previous case. In both cases, the values 
of headway in the jam and free driving region are just the same as analytic results. 
The total number of cars in jam denoted by Nj is given by the following formula, 
which does not depend on the initial conditions: 

(N - Nj)ax f + NjAXj ~ L . (32) 

The numerical result of Nj is just the same as that of the analytic prediction. 

Fig.^| is the plot of car points in the phase space accumulated over the steps 
after the jam cluster becomes stable. The data points are just on the hysteresis 
loop, which is derived analytically. We note the hysteresis loop is just the same 
for each cluster in both simulations. Actually, this profile is very universal, which 
does not depend on initial conditions nor boundary conditions. This fact is easily 
understood by the derivation of analytic results. In conclusion, all simulation 
data shows that the jam flow is realized just as the analytic solution prospects. 

We should discuss about the stability of the jam flow solution as well as the 
homogeneous flow solutions, which can be considered as "free driving flow" with 
no jam. The stable cases of the homogeneous flow are trivial. If all cars stop with 
the headway less than d, this situation is stable. While, if all cars are moving 
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Figure 5: The simulation data of the distribution of headway of each car. 
The data agrees with the analytic results: Axj ~ 0.41 and Axf ~ 3.59. 



with the same velocity v max and with the headway larger than d, this flow is also 
stable. In contrast to the above two cases, if the headway of a car is more than 
d while the others are less than d, or the opposite case, such flows are unstable 
and the jam flow is organized with time development. And, we can observe the 
jam flow is stable, as an example in Fig.|]. 

To summarize, we have proposed the simplest model for queue dynamics, 
which has the travelling cluster solution (jam flow) derived analytically. In the 
solution, cars are moving just the same manner with the induced time delay. 
The amplitude of the cluster is characterized by this time delay, which is related 
to the response sensitivity. The cluster has the profile of a limit cycle dynamics. 
This shows the "delay of changing motion" of each car, which means the balance 
of car moving in and out of a cluster, and this microscopic balance results the 
formation of the macroscopic cluster. Our simple model can provide the general 
understanding for queue dynamics by comparing with another approaches ||. 

We thanks K. Nagel and M. Schreckenberg for inspiring us the idea of this 
work through the discussion during the one of the author (Y.S) was staying in 
Duisburg. And we should comment that T. Nagatani has the work concerned 
with the similar simplest cellular automaton model based on the OV model . 
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